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Abstract. In this paper, so-called alternative mass balance equations for 

atmospheric constituents published recently are assessed in comparison with the true 
local mass balance equations deduced from exact integral formulations. It is shown that 
these alternative expressions appreciably violate the physical law of the conservation 
of mass as expressed by the equation of continuity. It is also shown that terms of these 
alternative mass balance equations have different physical units, a clear indication that 
these alternative expressions are incorrect. Furthermore, it is argued that in the case 
of alternative mass balance equations a real basis for Monin-Obukhov similarity laws 
does not exist. These similarity laws are customarily used to determine the turbulent 
fluxes of momentum, sensible heat and matter in the so-called atmospheric surface layer 
over even terrain. Moreover, based on exact integral formulations a globally averaged 
mass balance equation for trace species is derived. It is applied to discuss the budget 
of carbon dioxide on the basis of the globally averaged natural and anthropogenic 
emissions and the globally averaged uptake caused by the terrestrial biosphere and the 
oceans. 
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1. Introduction 

Recently, Finnigan et al. pjj and many others [21[Ill[I51[I71Enilinill21llllE2J proposed 
alternative forms of balance equations of atmospheric constituents. Unfortunately, these 
equations are not self-consistent and their terms are afflicted with different physical 
units. The latter is a clear indication that something is going wrong. Since mainly the 
results of turbulent fluxes of carbon dioxide (CO2) published recently were harmfully 
affected by these alternative mass balance equations, these results must be assessed 
carefully, especially in front of the debate on global warming due to an increase of 
greenhouse gases like CO2. 

To assess these alternative mass balance equations it is indispensable to compare 
them with the correct ones. Therefore, in section 2 we will derive the correct local 
mass balance equations from exact integral mass balance equations and discuss them 
physically. It is shown in section 3 that one of the fundamental physical principles 
namely the conservation of mass (as expressed by the equation of continuity) is notably 
violated by the alternative mass balance equations. Some aspects of Reynolds averaging 
are briefly, but thoroughly discussed in section 4. In identifying an adequate averaging 
time, we argue that the Allan variance criterion [1] combined with the wavelet analysis is 
much more favorable than block averaging as suggested, for instance, by Finnigan et al. 
[15] and Trevio and Andreas [M]- In section 5, we present the Reynolds averaged mass 
balance equations for a turbulent fluid and compare it with that of the alternative forms. 
In section 6, these different forms of mass balance equations are vertically integrated 
by assuming horizontal homogeneity and the results obtained are assessed. We show 
in section 7 that in the case of the alternative form of the mass balance equations 
a real basis for Monin-Obukhov similarity laws does not exist. Consequently, such 
similarity laws customarily used for determining the flux densities (hereafter, a flux 
density is simply called a flux) of long-lived trace gases like CO2 would be obsolete. 
Based on exact integral formulations a globally averaged mass balance equation for 
trace species is derived in section 8 to underline that the emission and uptake of matter 
serve as boundary conditions. This result is in complete contrast to the alternative 
mass balance equation in which emission and uptake are considered as volume-related 
source/sink functions. Final remarks and our conclusions are listed in section 9. 

2. The mass balance equations of atmospheric constituents for a 
macroscopic fluid 

First, we will derive the conservation laws for a macroscopic fluid (often called a 
molecular flow). The total rate of change in the mass Mj of a constituent, i = 0, 1, . . . , A^, 
within an arbitrary volume V{t) that also depends on time, t, can be expressed by [10] 




(1) 
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Figure 1. Transfer of matter through the surface A{V{t)) of a given volume V{t). 

where pi = Mi/V is the corresponding partial density. This quantity dMi/dt is equal to 
the material flow of the i*'' component into and out of this volume through its surface 
A{y{t)), i.e., the exchange of material between the volume and its surroundings, and 
the net production/depletion ai = /(pi, . . . , Pn) of the i^^ constituent per unit volume 
owing to chemical reactions (and phase transition effects) that occur inside V. Thus, 
we have (see Figure [T]) 

^ [ p,dV = - [ Pi ■ ciA + / a,dV . (2) 

Jv{t) JA(V{t)) JV{t) 

Here, v^v = Vj— v is the net velocity that characterizes the diffusion of the i^^ constituent 
with respect to the barycentric flow, considered as being of Newtonian kind. The 
quantity Vj is the individual velocity, and v denotes the barycentric velocity defined 
by 

1 ^ 

^ = - ^Pt^i , (3) 



=0 



where 



M 

P 



V V 



^ N N 

r,Y.^^ = Y.p^ (4) 

is the total density of air. The use of the barycentric velocity as a reference velocity, 
of course, is not the only possibility to describe diffusion. Prigogine [57] deduced from 
the entropy principle that in systems with mechanical equilibrium diffusion processes 
can be related to an arbitrary reference velocity. Herbert [231 El] discussed the 
general application of Prigogine's diffusion theorem to the atmosphere and some specific 
invariance properties of the thermodynamic laws as well as various alternative relations 
to describe Fick-type mass diffusion in a (diluted) binary gas mixture such as the 
atmosphere. Nevertheless, in our contribution the diffusion is related to the barycentric 
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velocity and the quantity Jj = piVN = Pi{^i — v) is denoted as the diffusion ffux of 
the i^^ constituent. Furthermore, dA = n dA is a vector normal to the surface of the 
volume V with the unit vector n and the magnitude dA. The unit vector is counted 
positive from inside to outside of the volume. Note that i = denotes that portion 
of dry air which is chemically inert, i.e., this portion does not contain any chemically 
active atmospheric constituent. Furthermore, the occurrence of the surface integral in 
Eq. ([2]) means that we consider a system that is open in the sense of thermodynamics, 
when the exchange of energy is allowed, too. It is the most general one. 

Since the volume is considered as time-dependent, the differentiation of the volume 
with respect to time is also required. It can be performed by applying Leibnitz's integral 
theorem. In doing so, we obtain for the total rate of change, dMi/dt, 

tJ P^dV= [ ^dV+[ p.vdA. (5) 

(^t Jv(t) JV(t) JA(V(t)) 

Thus, combining equations ^ and ([S]) results in 

JV(t) Jv{t) JA{V{t)) 

Ji-dA+ [ aidV (6) 

AiVit)) Jv{t) 



or 



/ 7^^^+/ (piV + J,)-t/A= / a,dV. (7) 

Jv{t) JA(V(t)) Jv{t) 

This expression is called the integral balance equation of the i^^ atmospheric constituent. 
It is a very general formulation and not affected by any kind of averaging process usually 
applied in deriving local balance equations for turbulent systems like the atmospheric 
boundary layer. McRae and Russell |15], for instance, used such an integral formulation 
for practical purposes by considering the domain of the South Coast Air Basin of 
Southern California. 

By applying Gauss' integral theorem (often called the divergence theorem), the 
surface integral occurring in Eq. ([T]) can by replaced by a volume integral so that we 
obtain 

'^dV+ [ V ■{piY + 3,)dV = [ aidV, (8) 

V{t) Ot Jy(^^-j Jy(^f.^ 



where V is the nabla (or del) operator. Since in Eq. ([8]) the shape of the volume is 
arbitrary, we may also write: 

^ + V-(AV + J,) = cr,. (9) 

This equation is called the local balance equation of the i*'^ constituent P, [101 ESI EH EH 
[M] . The term dpi/dt is the local temporal change of pi, and the quantities pjV and Jj are 
frequently called the convective and non-convective transports of matter, respectively. It 
is obvious that in the case of a closed or isolated thermodynamic system the divergence 
term vanishes because it only describes the exchange of matter between the system and 
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its surroundings. As the use of Gauss' integral theorem requires that some mathematical 
pre-requisites have to be fulfilled, the local balance equation ^ is somewhat lesser valid 
than the integral balance equation ([t]). Nevertheless, any integration of Eq. ^ over a 
certain volume as performed, for instance, by numerical models of the troposphere must 
be in agreement with Eq. ([T]). 

Summing Eq. ^ over all substances, i = 0, 1. . . . , A^, provides 

N N N 

E|^+Ev-(p.v+j.)=j:a.. (10) 

i=0 i=0 i=0 

Since, according to de Lavoisier's law, mass is conserved in chemical reactions (and/or 
phase transition processes; i.e., the transformation of mass into energy as expressed by 
Einstein's formula is unimportant in the case of atmospheric trace constituents), we 
have 

N 

^a, = 0. (11) 

i=0 

Furthermore, with the aid of Eqs. (|3]) and ^ and the definition of the diffusion flux, 
we obtain 

Et^|(E.M (-) 

1=0 i=0 



and 



N N 



^ V • (pi V + J,) = V ■ (^(p, V + p,(v, - v))^ 

i=0 i=0 

N 

= V-(5^(p.v,))=V-(pv). (13) 

2 = 

Thus, Eq, (10) leads to the macroscopic balance equation for the total mass per unit 

volume [ini [201 E21 EZl EHl ED 

|^ + V-(pv) = 0. (14) 

It represents the formulation of the conservation of mass on the local scale and is 
customarily called the equation of continuity. 

If we use the mass fraction q = Pi/p, then Eq. ^ will read 

^^ + V-(pc,v + Ji) = a,. (15) 
Applying the product rule of differentiation to this equation yields 

P ^ + Ci(|^ + V • (p v)) + p V ■ Vq + V ■ Ji = a, (16) 

= (seeEq.^) 

or 

fdc, 



p(^^+v-Vcij +V-J, = ai , (17) 
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where 



dci 
~dt 



dcj 

at 



+ V ■ Vci 



is the substantial or total derivative with respect to time. Equation (17) may be called 
the advection- diffusion equation of a chemically active atmospheric constituent. 



3. The alternative mass balance equations of atmospheric constituents 

Recently, several authors proposed a mass balance as an alternative form to the 
conservation equation (g . This alternative form reads [2l[131[T5|[T6|[T71|89l[inl[i2|[lit|62] 

^ + V-(p.v) = 5.. (19) 

In this equation, the diffusion flux Jj is generally ignored. This flux, however, is 
not negligible because, as expressed in Eq. (|9]), it must represent the local balance 
equation for a macroscopic fluid. If Jj would not occur, diffusion of the z*'* constituent 
in a macroscopic fluid could not be quantified. This means, for instance, that the 
sedimentation of airborne particles like aerosol and ice particles and water drops would 
generally be excluded. Moreover, this alternative form is applied to quantify the sinks 
or sources of CO2 inside canopies of tall vegetation like forests by a method where the 
term Si is called a source/sink inside the volume V of the fluid. However, considering 
the derivation of Eq. (|9]), an uptake or emission of a substance by plant elements or soil 
is a process at the boundary of a fluid and, therefore, has to be described in terms of 
boundary conditions, as substantiated by the integral balance equation ([t]). Therefore, 
Si is not a biological source/sink term. Such a term must not occur in the local form 
of a balance equation. Boundary conditions only occur when the local mass balance 
equation is integrated, in complete agreement with Eq. ([T]). The alternative equation 



(19), however, cannot be deduced from any integral mass balance equation. Finnigan et 
al. jl5] introduced it into the literature in an unforced manner and without any physical 
justification. One of the physical consequences related to this alternative mass balance 
equation is that the biological source/sink would explicitly cause a temporal change in 
the partial density. Whereas Eq. ^ clearly substantiates that only the divergence of 
the convective and non-convective fluxes contributes to a temporal change of the partial 
density. 

Finnigan et al. [15] and many others |2l[Tlll6l|l71|39llinill2lll2[62] expressed Si 
as a surface source (instead of a volume-related source or sink due to chemical reactions 
or /and phase transition processes) by 

5, = 5,(x)5(x-xo), (20) 

where 5(x — xq) is Dirac's delta function (from a mathematical point of view a 
distribution). Finnigan et al. [15] argued that the source term, (x), is multiplied 
by the Dirac delta function, signifying that the source is zero except on the ground and 
vegetation surfaces, whose locus is xq. This is incorrect at least by two reasons: As 
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mentioned before, the physical meaning of 5*^ in Eq. ( 19 ) must be a source or sink within 
the volume, but not at the boundaries, which are represented by boundary conditions. 
Furthermore, the Dirac delta function has the property that 5(x — xq) = for all x 7^ xq 
[5l [m [381 SSI EH], i.e., there is only the infinitesimal region (namely when x = xq) at 
which a source or sink of matter is defined. In addition, in various textbooks [SSI SS] , 
it is pointed out that (5(x — xq) = 00 for x = xq. Thus, Eq. (20) and, hence, Eq. (19) 
would become meaningless for x = xq. As pointed out later on, the main properties of 
the (5-function are defined by the integral over a region containing x = xq. 

Uptake and emission of a constituent by plants, soil and/or water systems are 
surface effects, expressed, for instance, in SI units by kg/{rn?s)] whereas the local 
temporal change in Eqs. ^ and (15) has the SI units kg/{m^s). If we express Dirac's 
delta function using SI units we will have m~^. Therefore, the term ^^(x) 5(x — xq) 
would have the SI units kg/{m^s). It is indispensable that in any physical equation its 
terms must have identical physical units. This fundamental requirement, however, is 
not satisfied in the case of the alternative forms of mass balance equations. There is 
evidence (e.g., Eq. (6) of Finnigan et al. [15], Eq. (11) of Finnigan [16], Eqs. (2.2) and 
(3.8) of Lee et al. [10], and Eq. (14) of Massman and Tuovinen [H]) that the biological 
source/sink term is, indeed, expressed by the units of a flux of matter, i.e., in kg/ {m?s), 
even though, as expressed by the local temporal change of the partial density (dpi/dt), 
kg/{m?s) is unequivocally required. 

Nevertheless, assuming for a moment that Eq. (20) is reasonable, then Eq. (11) 
would result in 

N N 

^5. = ^^.(x)5(x-xo)^0, (21) 

i=l i=l 



and, hence, Eq. (10) in 

^ + V ■ (p v) = 5^ 5.(x) 5(x - xo) ^ . (22) 

i=l 

The right-hand side of this equation is not equal to zero because atmospheric 
constituents have different natural and anthropogenic origins and different sinks. This 
means: if Eq. (20) would be reasonable, than one of the fundamental laws of physics, 
namely the conservation of the total mass on a local scale, as expressed by the equation 
of continuity (14), is notably violated by the alternative forms of mass balance equations. 



4. Reynolds averaging calculus 

Conservation laws for a turbulent fluid can be derived by using Reynolds' averaging 
calculus, i.e., decomposition of any instantaneous field quantity ^{v) like Pi{Tc) and v(r) 
hj ip{v) =Tp + Lf' and subsequent averaging according to [22l [321 169] 

V = W) = ^lj{r,r')dG' , (23) 
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where Ip is the average in space and time of V5(r), and the fluctuation ip' is the difference 
between the former and the latter. Here, r is the four-dimensional vector of space and 
time in the original coordinate system, r' is that of the averaging domain G where its 
origin, r' = 0, is assumed to be r, and dG' = d^r dt'. The averaging domain G is given 
by G = J^dG'. Hence, the quantity ^ represents the mean values for the averaging 
domain G at the location r. Since Tp = Tp^ averaging the quantity ipij) =ip + ({)' provides 
ip' = 0. Any kind of averaging must be in agreement with these basic definitions. These 
basic definitions cannot be undermined by imperfect averaging procedures. Since, in 
accord with the ergodic theorem, assemble averaging as expressed by Eq. (23) may be 
replaced by time averaging [H], the common practice, it is indispensable to ask whether 
a time averaging procedure is in complete agreement with the statistical description of 
turbulence or not. 

A basic requirement for using a time averaging procedure is that turbulence is 
statistically steady [131 HH [67] . Consequently, sophisticated procedures for identifying 
non-stationary effects (trends) are indispensable to prevent that computed turbulent 
fluxes of atmospheric constituents are notably affected by non-stationary effects. In 
order to identify such stationary states in the off-line time series analysis, Werle et al. 
[70] suggested the Allan-variance criterion [1] as a suitable and efficient tool. As argued 
by Percival and Guttorp [31] , this variance is an appropriate measure for studying long- 
memory processes because it can be estimated without bias and with good efficiency 
for such processes, and it may be interpreted as the Haar wavelet coefficient variance 
[T8l l3T] . Percival and Guttorp [5l] generalized this to other wavelets [36] • Following 
various authors [6l EU [261 12H] , wavelet decomposition seems to be well appropriate to 
study turbulence data. Thus, the Allan variance criterion combined with the wavelet 
analysis is much more favorable than simple block averaging as suggested, for instance, 
by Finnigan et al. [15] and Trevio and Andreas |68j . 

Figure [2] illustrates two Allan plots deduced from a synthetic data set generated 
by the equation Yt = a + b t + G{a) (see Figure [3]), where a is an offset, 6 t is a 
linear drift and G{a) is a Gaussian distributed white noise (lower trace) and a data set 
containing a ten times stronger drift (upper trace). In both cases CF\{k) was calculated 
using Haar wavelet coefficients [31j. For convenience, it is assumed that the sampling 
interval between consecutive observations is constant and amounts to At = 1 s. Thus, 
the time r = /c At is equivalent to the averaging or integration time. As shown in 
Figure [3| the lower trace indicates a minimum of (j\{k) at an integration time, denoted 
here as optimum integration time Topt-, of about 130 s. Whereas the upper trace for 
the data set with the ten times stronger drift suggests a Topt value of about 25 s. For 
T < 'Topti the Allan plots show a behavior that is typical for white noise. Beyond Topt-, 
the Allan plots obey the law which indicates that a linear drift becomes dominant. 
Consequently, for r > Topt, stationary conditions as required by time averaging are not 
assured. From this point of view, Topt may be considered as the maximum averaging 
time 
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Figure 2. Allan plots deduced from the time series data illustrated in Figure |3] (solid 
line) and a data set which differs from that by a ten times stronger drift (dotted line). 
The dashed lines indicate white noise (r^^) and linear drift (r^) behavior (adopted 
from f31j). 



5. The balance equations of atmospheric constituents for a turbulent fluid 

Applying the conventional Reynolds' averaging calculus to Eq. ^ yields [HI 



+ V- (^p7v + p', v' + J,j =(T^. (24) 
Here, the overbar represents the conventional Reynolds average; whereas the prime 



denotes the departure from that. Equation (24) is called the balance equation for the 
1** moment (or l**-order balance equation). Here, p[ V is the turbulent (or eddy) flux. 
It is a non-convective flux, too. 

As mentioned before, the diffusion flux represents not only molecular fluxes, but 
also the sedimentation of aerosol particles affected by the gravity fleld. Thus, Jj cannot 
generally be neglected in comparison with the corresponding second moment p[ v'. Even 
in the immediate vicinity of the earth's surface it must not be ignored because the wind 
vector vanishes at any rigid surface, and, consequently, no exhalation or deposition of 
matter would be possible. In the case of gaseous constituents Jj becomes negligible in 



comparison with v' when a fully turbulent flow is considered. Obviously, Eq. (24) 



also fulflls the conditions Yld=Q Jj = and Yl!i=i — 0' ^^"^^ required by the equation 
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Figure 3. Generated time series data of a concentration consisting of an offset, a 
linear drift and a Gaussian distributed white noise. The sample size amounts to 
iV = 2i" = 1024 (adopted from [7T]). 



of continuity in its averaged form, 

In 



+ V ■ p V + p' 



(25) 



the conditions Yld=o Pi v = p v and ^ P'i ^ 

If mass fractions are considered (see Eqs. (15) to ([li 
averaging techniques [121 1221 EOl ESI EH] should be apphed in formulating local balance 
equations for turbulent flows to avoid any simplification |32l [3] 



N 



p' v' are fulfilled, too. 

then density-weighted 



If we rearrange Eq. (24) by using the equation of continuity (25) (it is similar 
to the rearrangement of Eq. (15) that leads to Eq. (17)) and ignoring the covariance 
terms p' v' and p' c'^, as customarily accepted within the framework of the Boussinesq 
approximation, we will obtain 
.(del 



V • 



V 



(26) 



/^ V + 

It may be called the advection-diffusion equation of a chemically active atmospheric 
constituent for a turbulent fluid. Here, cTj still represents the gain or lost of matter due 
to chemical reactions (and/or phase transition processes). 

For the purpose of comparison: Finnigan's [TB] Eqs. (8) as well as Eq. (3.7) and 
(10.8) of Lee et al. [40] read 



dt 



+ V ■ Vp7 + V ■ p 



(27) 
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where Si is still given by S'i(x)5(x — xq) (see our Eq. (20)). Obviously, Eq. (27) is 
not self-consistent. First, as shown in Eq. (26), the mass fraction Cj must occur in the 
two terms on the left-hand side of Eq. (27), but not the partial density pi. Second, the 
rearrangement of Eq. (15) that leads to Eq. (17) requires that the right-hand side of 
Eq. (22) is always equal to zero so that Eq. (22) equals the equation of continuity (14). 



6. Vertically integrated mass balance equations 



If we assume horizontally homogeneous conditions and recognize that in such a case the 



vertical component of the mean wind vector becomes nearly equal to zero 



w 



(24) will provide 
dt 



+ 



d_ 
dz 



p[ w' + Ji,, 



0), Eq. 



(28) 



where, w is the vertical component of the wind vector, Jj^^ is the vertical component of 
the diffusion flux, and z is the height above ground. Note that Archimedean effects are 
related to the gravity field. The same is true in the case of the hydrostatic pressure. 
Consequently, the vertical direction is related to the gravity field[nS]; but not to the 
normal vector of a slope or a streamline as recently described, for instance, by Finnigan 
et al. [IS] and Finnigan (TU]. Especially over complex terrain trajectories and streamlines 
are usually vary with time and position. Therefore, it is difficult to find, for instance, a 
common reference streamline coordinate frame for using sonic anemometry at different 
heights to estimate the variation of turbulent fluxes of momentum, sensible heat and 
matter with height. Without knowing a reference coordinate frame that is invariant 
in space and time, the convergence/divergence of any trace gases cannot be calculated 
[651 [66] . In contrast to a streamline coordinate frame, one may use a terrain-following 
coordinate frame as usually considered within the framework of mesoscale meteorological 
modeling, but an exact transformation of the governing equations is still indispensable. 
Fortunately, this transformation is well known since more than three decades (see, e.g., 

m)- 

The integration of this equation from the earth's surface {z = 0) to a certain height 
above ground {z = /i), where a fully turbulent flow is established, yields then 

"""'dpi 



dt 



a, 



dz 



p\ w' + Ji^^ 



Pi W + Ji, 



Assuming long-lived trace species (dj = 0) and stationary condition {dpi/dt 
to the constant flux approximation expressed by 



p[ w' + Ji, 



P'i W + Ji, 



(29) 
0) leads 

(30) 



At the height z = h the vertical component of the diffusion flux of a trace gas can usually 
be ignored in comparison with the vertical component of the corresponding eddy flux 
component {\Ji^z\z=h |p'j w'\z=h). As already mentioned, the opposite is true in the 
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immediate vicinity of the earth's surface {\Ji^z\z=o ^ \p'i w'\z=q). Thus, Eqs. (29) and 

a-i]dz= J,., - p',w' (31) 



(30) may be written as 



dt 



2=0 



and 



p'i w 



Ji. 



J 2=0 



(32) 



In contrast to this, Eq. ^ with Si = Si{z) 6{z) (see Eq. ([20])) leads to [ISl [161 El HT] 

(33) 



Here, an important inconsistency exists because x and xq are vectors so that (5(x — xq) = 
has to be expressed, for instance, in Cartesian coordinates by 6{x—xo)6{y—yo)6{z) = 
when = is considered. Therefore, we would have 



+ §^{p'i = ^i(^) ^o) S{y- yo) 6{z) 



(34) 



Nevertheless, following for a moment Finnigan and disciples, the integration of this 
equation should yield [T^ [TBI [TTlliT] 

"''dpi 



dt 



-dz 



p'i w 



where the biological source/sink term Si is given by 

_ 

Si = Si{z) 6{z) dz . 
Jo 

Assuming steady-state condition yields then |17] 



p'i w 



Si . 



(35) 



(36) 



(37) 



Obviously, Eq (37) looks similar like Eq. (32). However, from a mathematical point of 



view, Eq. ([36[) is faulty because Dirac's delta function has the fundamental property 

(38) 



that [in [Miiini [ss 

f{x) 5{x — a) dx = f{a) , 
where f{x) is any function continuous at the point x = a, and, in fact. 



f{x) 6{x — a) dx = f{a) for e > 



or, with bi = a — e < a and b2 = a + e > a, 

/ f{x) 6{x — a) dx = f{a) for 6i < a < 62 ■ 
In the special case of /(x) = 1 this equation directly provides 
5{x — a) dx = 1 . 



(39) 



(40) 



(41) 
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This means that the range of integration must include the point expressed in 



Eq. (40) by 61 < a < 62; otherwise, the integral equals zero. Equation (36) does not 
fulfill this requirement. Consequently, we have 



S, = / Si{z) 6{z) dz = . (42) 
Jo 

This result is generally valid for any value of z > 0. 



7. The bases for Monin-Obukhov similarity laws 

Similarity hypotheses according to Monin and Obukhov [3l HH] are based on the pre- 
requisite that the turbulent fiuxes of momentum, sensible heat and matter are invariant 
with height across the atmospheric surface layer (ASL). These similarity hypotheses can 
be expressed by [Ml El [521 129] 



and 



~ = H0- (45) 

C-i,* QjZ 



Here, k is the von Krmn constant, u^, is the friction velocity, U = |v/7| is the magnitude 
of the mean horizontal wind component, is the potential temperature, 0* is the so- 
called temperature scale, and q,,, is the so-called scale of matter. Furthermore, $m(C); 
$/i(C), and $j(C) are the local similarity functions for momentum (subscript m), sensible 
heat (subscript h), and matter (subscript i), respectively. Moreover, ( is the Obukhov 
number, where L is the Obukhov stability length defined by [Sni SHI [121 [22] 

L = *- = . (46) 

I (0* + 0.61 c,,,) 

Here, g is the acceleration of gravity, and c^^^: is the so-called humidity scale when 
subscript i stands for water vapor. The turbulent fiuxes of momentum (the magnitude 
of the Reynolds stress vector), r, sensible heat, H, and matter, Fj, are related to these 
scahng quantities bv pi [8l [521129] 



T = [ p u' w' + p v' w' ) = p = const. , (47) 



H = Cp^cL p w' &' = —Cp d p = const. , (4J 



and 



Fi = p w' c'j = —p Cj_=K = const. , (49) 
where Cp^d is the specific heat at constant pressure for dry air. Assuming that these 



fiuxes are invariant with height and integrating Eqs. (43) to (45) over he height interval 
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[zi, Z2\ {zi and Z2 may be the lower and upper boundaries of the fully turbulent part of 
the ASL) result in 



and 
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Ciizi) 
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'in 


K ' 






(in 


K 






(in 







--^.(C2,Ci; 
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Zl 



where the corresponding integral similarity functions are defined by [5T 



■m,h,i 



(C2,Cl 



Ci 



C 



dc 



(50) 
(51) 

(52) 
(53) 



Obviously, not only the similarity hypotheses of Monin and Obukhov but also the 
integration of Eqs. (43) to (45 require that the fluxes of momentum, sensible heat 
and matter are height-invariant across the ASL. Since the local similarity functions 
of Monin and Obukhov cannot be derived using methods of dimensional analysis 
like Buckingham's [7j vr theorem, it is indispensable to use formulae empirically 
derived. Integral similarity functions that are based on various local similarity functions 
empirically determined are reviewed by Panofsky and Button [32], Kramm [33], and 
Kramm and Herbert l35l. 



The set of so-called profile functions (50) to (53) enables the experimentalist to 
estimate turbulent fluxes of momentum, sensible heat and matter in dependence on 
the thermal stratification of air at a certain height above the lower boundary, namely 
the earth's surface, using the vertical profiles of mean values of windspeed, potential 
temperature and mass fractions [201 Furthermore, the parameterization schemes 
used in state-of-the-art numerical models of the atmosphere (weather prediction models, 
climate models, etc.) for predicting the exchange of momentum, sensible heat and 
matter between the atmosphere and the underlying ground also based on this set of 
profile functions. 

In the case of trace species Eq. (24) is the essential rule for micrometeorological 
profile measurement. With respect to this, it requires stationary state and horizontally 
homogeneous conditions and that chemical reactions play no role like in the case of CO2. 
Under these premises one obtains 
d 



dz 



Pi W + J^, 



^ p'j w' + Ji^z = const. 



(54) 



This result completely agrees with Eq. (30). On the contrary, under the same premises 
the alternative equation (27) leads to 
d 



dz 



Pi w 



Si 6{z) 7^ const. 



when we assume for a moment that Eq. (36) would be correct. 



(55) 

Consequently, the 



alternative advection-diffusion equation (27) would imply that no basis for the similarity 
laws of Monin and Obukhov does exist, i.e., the profile functions (50) to (53) customarily 
used for determining the fluxes of long-lived trace gases like CO2 would be obsolete. 
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Figure 4. Atmospheric COi concentration inferred from Siple station (West 
Antarctica) ice cores analyzed by Neftel et al. [49 and Friedii et al. [19 and extended 
to that directly observed at various locations during the period 1958-2004 (courtesy 
of CD. Keeling, T.P. Whorf, and the Carbon Dioxide Research Group of the Scripps 
Institution of Oceanography (SIO). The fit is based of a polynomial of S*'' order. 



8. The global budget of carbon dioxide 

Dividing Eq. ([T]) by tlie volume and rearranging the resulting equation yield for a 
turbulent system 

'%)=-ViT\ [ i7h.^ + P^ + Jl)-dA+{al)y, (56) 
at / V V [t) jAivm 



my = - I ijV. (57) 



'A{V{t)) 

where the volume average of an arbitrary quantity ip is generally defined by 

1 

V 

Hitherto, the volume V(t) has been considered as arbitrary. Suppose the earth can be 
considered as a sphere with the radius te = 6, 371 km and the atmospheric layer under 
study that directly covers the earth as a spherical shell of a thickness of Ar, the volume 
of the atmospheric layer, Va, now considered as independent of time is then given by 

Since Ar/rE ^ 1, the term (1 + At/teY can be approximated by 1 + 3 Ar /te- Thus, 
we obtain 

= 4 vr r| Ar . (59) 
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The surface of this spherical shell is given by Aa = At + Ae = 4 it r"^ {(1 + Ar/r£;)^ + 2}, 
where At and Ae are the outer surface and the inner surface of this shell, respectively. 



The latter is congruent with the earth's surface. The surface integral in Eq. (56) may, 
therefore, be expressed by 



{-pi^ + Pi v' + Ji) -rfA 



Aa 



{pi V + p'i v' + Ji) ■ dA 



At 



+ / {p-iV + p'-V + Ji)-dA 



(60) 



The first integral on the right-hand side of this equation describes the exchange of 
matter between the atmospheric layer under study and the atmospheric region aloft 
across the common border, and the second integral the exchange of matter between 
the atmospheric layer under study and the earth caused by the total (anthropogenic 
plus natural) emission (counted positive) and the uptake (counted negative) by the 
terrestrial biosphere (plants and soils) and the oceans. If we choose Ar in such a sense 
that exchange across the outer surface of this spherical shell is equal to zero (or in 
comparison with that at the earth's surface, at least, negligible) because of the vertical 
profile of the CO2 concentration becomes (nearly) independent of height (see, e.g.. 
Figure 1 in [63]), we will obtain 



{p-iY + p'- V + Ji)-dA= / (p- V + p'. v' + Ji) ■ c/A . (61) 

'Aa J Ae 

Since the unit vector normal to the inner surface of the spherical shell shows in the 
direction of the earth's center, we have only to consider the radial components of these 
fluxes characterized by the subscript r, i.e.. 



{pi V + p'. V' + Ji ) ■ C?A = - / (pi Vr + p'iV'^ + Jr,i ) 
Ae J Ae 



dA 



EM 



+ / {pi Vr + PiK + Jr,i ) 
'Ae 



dA . 



U 



(62) 



The first integral on the right-hand side of this equation represents the total emission 
(subscript EM) and the second one the uptake (subscript U). The signs of these terms 
are determined by the scalar product between the unit vectors of the fluxes and the unit 
vector normal to the inner surface of the spherical shell. 



Combining Eqs. (56) and (60) to (62) yields 

/&P^\ ^ J_ 
\ dt Iva Va 



{Pi Vr + p'i V'^ + Jr,i ) 



dA 



EM 



{pi Vr + p'i < + Jr,j ) dA \ + {(Ti)v. ■ 
Ae ^ ^ 



(63) 



Since CO2 is a long-lived trace gas, the effect caused by chemical reactions can be 



ignored. Thus, Eq. (63) can be approximated by 



idp^ 

\ dt I Va 



1 

Va 



{pi Vr + Pi V'^ + Jr^i 



Af 



dA 



EM 
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dA 



E 



(Pi Vr + p'i < + Jr,i ) 

U J 

(64) 

Expanding the right-hand side of this equation with Ae yields finally 



EMI Ae 



Ul A 



{Pi Vr + p- < + Jr,i ) ) L (65) 



where 



{i^)AE = ^ I ^dA. (66) 



E J A 



represents the global average of an arbitrary quantity ip. Equation (65 ) deduced by using 
the exact integral formulations substantiates that (a) total emission and the uptake of 
CO2 have to be considered as lower boundary conditions and (b) the partial density of 
CO2 averaged over the volume of the atmospheric layer under study will rise as long as 
the globally averaged total emission is higher than the globally averaged uptake. 

9. Final remarks and conclusions 

In 1999 Sarmiento and Gruber [60] pointed out that "the land sink for carbon is 
the subject of considerable controversy at present, concerning not only its magnitude 
but also its cause". It seems that any use of an alternative mass balance equation 
in micrometeorology may contribute to considerably more confusion because this 
alternative expression is clearly incorrect. 

The use of an alternative mass balance equation can harmfully affect not only the 
whole atmospheric budget of CO2, but also that of other greenhouse gases like water 
vapor (the most important greenhouse gas), nitrous oxide {N2O), methane {CH4), and 
ozone (O3). 



If Eqs. (19) and (27) are correct, the biological source/sink, for instance, would 
explicitly cause a temporal change in the partial density. However, the reality is different. 
Figure |4] illustrates that the atmospheric CO2 concentration has been rising since the 
beginning of the 18th century. If these results are correct as stated in various reports 



of the Intergovernmental Panel on Climate Change (IPCC), we can infer from Eq. (65) 
that the total emission of CO2 has always been higher than the CO2 uptake during 
the period covered by these results. Thus, lowering, for instance, the anthropogenic 
emissions of CO2 to those of the year 1990 would not reduce the CO2 concentration 
in the atmosphere (see Special Report on Emissions Scenarios (SRES) of the Working 
Group III of the Intergovernmental Panel on Climate Change (IPCC) and the IPCC 
IS92, |25ll2Z|)- However, if the CO2 uptake would rise due to a higher atmospheric CO2 
concentration, a stabilization of this concentration at a level of nearly 550 ppmV, i.e., 
higher than that of 1990, might be possible [25], l27] . 

A decrease of the atmospheric CO2 concentration can only be achieved when the 
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CO2 uptake by the terrestrial biosphere and the ocean is higher than the total emission 
CO2. As indicated by Beck's[l] inventory for the past 180 years that is based on 
more than 90,000 observations using chemical methods, there were CO2 concentrations 
appreciably higher than the current value of about 385 ppmV. Only when Beck's data 
are rehable, we may conclude that during various periods of the past the uptake was 
stronger than the total emission. 
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